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ABSTRACT 



A two-dimensional finite element model of a low frequency, 
free- flooded piezoelectric thin ring sonar transducer was 
developed for use with the ATILA code. 

Effective material properties for the ring were determined 
from in-air modal analyses using ATILA. These were adjusted to 
obtain the closest agreement between the calculated ring 
resonance frequency and coupling coefficient and measured data 
supplied by the manufacturer. In addition, equations for the 
electrical admittance for this transducer in air were 
developed based on an analytical approach published by Hong- 
zhang Wang (J. Acoust. Soc . Am. 79, 164-176 (1986)) . A Matlab 
program was written to implement these equations and to 
calculate various electromechanical network parameters. The 
calculated values agreed with measured values to within 5 
percent . 

An in-water harmonic analysis, in which the transmitting 
voltage response and the directivity pattern are computed, was 
performed using the ATILA code. In general, results obtained 
using the finite element model and manufacturer's measured 
values agree within 10 percent. The effects of changes in 
material properties on acoustic performance of the transducer 
were investigated; the results are discussed. 
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I. INTRODUCTION 



Free - flooded piezoelectric ring transducers have found 
increasing use in deep ocean sonar applications. Though simple 
in concept, the ability to reliably predict the performance of 
such devices is complicated by the nature of their 
construction. In particular, the acoustic performance of such 
transducers depends upon the physical properties of the 
encapsulant materials used, making their performance difficult 
to model . 

Standard multiport network analysis techniques, long used 
to model the electromechanical behavior of piezoelectric 
transducers, are usually restricted to one -dimensional motion. 
Newer types of transducer structures, such as the free -flooded 
ring, that have more than one important degree of freedom, do 
not lend themselves to this technique. Consequently, finite 
element analysis (FEA) techniques, which provide the 
capability to analyze engineering problems of many degrees of 
freedom, are frequently employed. 

The objectives of this research were (1) to numerically 
model a free -flooded piezoelectric ring sonar transducer and 
(2) using the developed model, predict the changes in acoustic 
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performance due to changes in the elastic properties of the 
encapsuiants used. 

The electroacoustic properties of this free- flooded 
piezoelectric ring sonar transducer were numerically modeled 
using the finite- element code "ATILA", developed by the 
Institut Superieur d' Electronique du Nord (Lille, France) 
[Ref. 1]. Since the ring is constructed of ceramic staves 
cemented together, the effective material properties for the 
ring were determined from in-air modal analyses using ATILA. 
The material property values were adjusted to obtain the 
closest agreement between the calculated ring assembly 
resonance frequency and coupling coefficient and measured data 
supplied by the manufacturer. 

In addition to the finite element analysis, equations for 
the electrical admittance of the thin ring transducer when 
unloaded were developed based on an analytical approach 
published by Hong- zhang Wang [Ref. 2] . A MATLAB program was 
written to implement and compute various electromechanical 
network parameters. The results obtained from the analytical 
method were compared with the measured values and those 
obtained by the finite element code. 

The finite element model built was further used to predict 
the changes in acoustical properties due to a change in the 
type of encapsulant materials used in manufacturing. 

The remainder of this thesis is divided into seven 
chapters. Chapter II describes the transducer under study. 



2 



Chapter III gives a general description of the finite element 
analysis and the standard network analysis of piezoelectric 
ring transducers. Chapter IV describes the development of the 
2 -dimensional finite element model of the transducer for use 
with the ATILA code. Chapter V presents and discusses the 
results of the in-air modal analyses and in-water harmonic 
analyses using ATILA, as well as the results of the analytical 
model. Chapter VI presents the conclusions. Chapter VII 
provides suggestions for future work. Appendix A contains a 
copy of the input data file for the finite element mesh used 
in the in-water harmonic analyses. The MATLAB program written 
to implement the analytical model for the transducer in air is 
attached as Appendix B. Appendix C provides a list of the 
typical values of material constants used in the modeling. 
Appendix D is a sample calculation of the effective elastic 
constants of a homogeneous ceramic ring incorporating the 
fiberglass wrap. 
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II. TRANSDUCER DESCRIPTION 



Figure 1 shows a schematic diagram of a free- flooded 
piezoelectric ring sonar transducer. When immersed in water 
and vibrating radially, such a transducer exhibits two 
fundamental resonant modes that can radiate acoustically into 
the water. They are the radial resonance of the ring and the 
cavity resonance of the enclosed water column. Radial motion 
of the cylinder walls will excite the symmetrical cavity modes 
in the enclosed water column. The radial resonance is the more 
important resonance. 




Figure 1: The Free Flooded Ring 
Transducer 
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The motor element (piezoelectric ring) of the transducer 
under study consists of a set of 72 tangentially polarized 
lead zirconate titanate, ?ZT, Navy Type I) ceramic staves as 
shown in Figure 2. 

PZT is an isotropic poivcrystaiine ceramic which exhibits 
piezoelectric properties after poling [Ref. 4] . It exhibits a 
high degree of symmetry in the plane transverse to the poling. 
These materials are usually treated as member of the 6 mm 
hexagonal class [Ref. 4] . 




Figure 2: The Ceramic Staves 



Each stave has dimensions (w x 1 x h) 0.25 x 0.6519 x 6.8 
in. (0.635 x 1.656 x 17.272 cm) . They are cemented together 
using epoxy resin to form a ring of 14.44 inches (0.3668 m) 
outside diameter, as shown in Figure 3. 
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The ceramic plates are poled tangentially and electrically 
wired in parallel. Thus, it is apparent that, during 
operation, the electric field will either simultaneouly be 
parallel or antiparallel to the poling direction shown in 
Figure 3 . The ring is wrapped on the outside by epoxy- 
impregnated fiberglass of 0.03 in (0.762 mm) thickness, which 
provides a compressive stress bias. 

A cutaway view is as shown in Figure 4 . 
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Figure 4: Cutaway View of Ring Transducer 

The piezoelectric ring and fiberglass wrap is enclosed by 
two encapsulant materials: the inner one is a polyurethane and 
the other is a neoprene. Of particular importance is the 
annular mounting plate at the bottom, which has an inner 
radius less than half the radius of the transducer. Because of 
the narrow size of the center hole of the mounting plate, it 
is expected that the mounting plate will influence the 
acoustic cavity mode of the enclosed water column in the 
middle of the transducer. 
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Appendix C lists the physical dimensions and constants for 
the materials of the transducer. 

The main operational characteristics of the transducer as 
provided by the manufacturer are: 

a. Resonance frequencies are 1.0 kHz (cavity) and 2.5 
kHz { ring ) ; 

b. Source Pressure Level, which corresponds to the 
effective pressure on the acoustic axis, at ring resonance 
is 132 dB ref 1 /i Pa/V at 1 meter; 

c. Operational Depth, which corresponds to the maximum 
depth where the performance of the transducer is not 
compromised - unlimited. 

The measured admittance circles of the ceramic ring 
assembly without and with fiberglass wrapping in air are as 
shown in Figures 5 and 6, respectively. 
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Figure 5: Measured Admittance Circle of the 
Unwrapped Ceramic Ring Assembly, In 
Air 
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Figure 6: Measured Admittance Circle of the 
Ring Assembly with Fiberglass 
Wrap, In Air 
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The measured transmitting voltage response (TVR) with 1 
V drive and the vertical directivity pattern at 2 kHz are as 



shown in Figures 7 and 8, respectively. 




At the ring resonance frequency of 2750 Hz, the TVR is 132 
dB; it is 127 dB at the cavity resonance frequency of 1250 Hz. 
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The beam pattern is omni-directional ^n the equitonal 
plane and toroidal in a plane passing through the axis of 
symmetry, which in Figure 8 below is the 90° and 270° line. The 
mounting plate is at 270°. The narrow lobe at 270° is probably 
due to the mounting plate. 




Figure 8: Vertical Directivity 
Pattern at 2 kHz 
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III. METHODS OF ANALYSIS 



A. INTRODUCTION 

This chapter presents the simple approximate expressions 
for the ring and cavity mode frequencies from G.W. McMahon 
[Ref. 3] and discusses the theoretical background behind the 
finite element code, ATILA, as well as the analytical 
approach, suggested by Hong -zhang Wang [Ref. 2] , to obtain the 
approximate electrical admittance equations for the 
tangentially polarized piezoelectric thin circular cylindrical 
tube in the unloaded case. 

B. RING AND CAVITY MODE FREQUENCIES [Ref. 3] 

The behavior of the cavity and ring resonance in small 
length- to- radius ratio open tubes of lead zirconate titanate 
ceramic immersed in water was studied by G.W. McMahon. 

When the tube is short i.e., h/a <4 n , G.W. McMahon gives 
the following formulas: 




( 1 ) 



r 



n 

pa* 



7TC 0 



( 2 ) 



c 



(h + 2 ota ) 
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where w r is che ring mode anaular frequency, w. is the angular 
frequency of the first cavity mode, a is the inside radius, h 
is the height of the ring, 3 is the bulk modulus of water, t 
is the wall thickness of ring and c is the speed of sound in 
water. c 0 is the velocity of sound in the water column given 
by Equation 5 below and a, an end correction, is given by the 
empirical formula: 

a = C .633 - 0 . 106 Q (3) 



The dimensionless frequency parameter, Q = o).a/c f)l is given by 



Q| — -0.633j-0.106 0 2 = 1 (4) 

\ 2 a / 2 



The value of 0 is between 0.33 and 3.3. The velocity of sound 
within the tube, c 0 , is effectively less than the velocity of 
sound c in open water because of the finite stiffness of the 
tube wall . It is given by 



c 0 = c 



1 + 



2 Ba 

Y33 1 



(5) 



The calculated value of the cavity mode frequency for the 
transducer is presented in Chapter V, Section C. 
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C. FINITE ELEMENT ANALYSIS, THE ATILA CODE 

The application of finite element analysis (FEA) to solve 
boundary value problems consists of the transformation of the 
governing differential or integral equation (s) into a multi- 
nodal matrix equation, the solution of which represents the 
discretized solution of the problem. There are many techniques 
to obtain a finite element formulation [Ref. 5, 6, 7, 8] . 

ATILA is a finite element code developed at l'Institut 
Superieur d' Electronique du Nord (ISEN) in France for the 
analysis of underwater transducers. It utilizes the 
variational formulation of the finite element problem [Ref. 
9,10,11,12,13] . 

ATILA uses quadratic isoparametric elements. Isoparametric 
means the same polynomial (quadratic in this case) is used to 
interpolate both the geometry and the field variation. ATILA 
has 46 different types of elements, including shell, plate, 
transition, spring, trilaminar, and two- and three-dimensional 
solid elements of various types. With ATILA, it is possible to 
model lossy elastic, piezoelectric, magnetostrictive , magnetic 
and composite materials, fluid, solid- fluid interfaces, 
dampers, and radiation conditions. 

ATILA can perform: (1) static analyses, (2) modal 
analyses, which correspond to a free vibration problem, where 
the eigenf requencies and eigenmodes are computed, and (3) 
harmonic analyses of radiation or scattering problems, which 
correspond to a forced vibration problem, the excitation being 
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che vo 1 cage applied across che electrical termmaxs of the 
transducer or external forces applied to the nodes. 

1. Harmonic Analysis of a Radiating Piezoelectric 

Transducer 

This problem is governed by the equation of motion in 
the elastic and piezoelectric structures, by Poisson's 
equation in the piezoelectric structures, and by Helmholtz's 
equation in the fluid. Appropriate boundary conditions are 
defined, both on the solid- fluid interface and on the external 
fluid boundary, which must simulate the appropriate radiation 
condition . 

The solid equation of motion is given by [Ref. 12] : 



P 



d 2 u f 

JF 






( 6 ) 



where p is the solid material density, u is the displacement 
vector, t is time, a is the stress tensor, and Xj is a 
coordinate direction. Here i and j can be 1, 2 and 3, and the 
Einstein notation is used, where summation is implied over 
repeated indices in the same term. 

Poisson's equation is given by [Ref. 12] : 



dD t 

bx t 



= 0 



(7) 
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where D is the electric displacement vector and is a 

coordinate direction; 1 can be 1 . 2 and 3. 

The linearized, lossless Helmholtz Equation for the 
propagation of sound in fluids is given by [Ref. 14] : 

V 2 p - 1 il = 0 (8) 

c~ d t" 



where V 2 is the three-dimensional Laplacian operator, p is the 
acoustic pressure, and t is time. 

In piezoelectric materials the stress tensor and the 
electric displacement vector can be related to the strain 
tensor and electric field vector via material properties by 
the following constitutive equations, which neglect magnetic 
and pyroelectric effects [Ref. 12]: 

T,J = c* m S u - e !Jk E k (9) 

D i = S a , Ej (10) 

where [T] is the stress tensor, [S] is the strain tensor, E is 
the electric field vector, D is the electric displacement 
vector, [c®] is the constant electric field elastic stiffness 
tensor, [e] is the piezoelectric tensor, and [ e 5 ] is the 
constant strain dielectric tensor; i, j, k and 1 can be equal 
to 1 , 2 and 3 . 
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Ultimately the solution must be expressed in terms of 
displacements and electric potentials. To this end the 
following two equations from elasticity and electricity, 
respectively, are used [Rer . 12] : 



S, 



^ i 

— \ 



<3 u 



du 

1H. 



(ID 






dx. 



( 12 ) 



where [S] is the strain tensor, u is the displacement vector, 
x± is a coordinate direction, E is the electric field vector 
and <f> is the electrical potential; i and j can be equal to 1, 
2 and 3. It must be noted that the quasi-static approximation 
[Ref. 4] is used in writing Equation 12 for the low (acoustic) 
frequencies involved here. 

The boundary conditions and prescribed excitations at 
each node can be defined either by a displacement or an 
applied force, an electrical potential or an electrical 
charge, or an acoustic pressure. 

In ATILA, the previous seven equations are transformed 



into the 


following 


matrix 


equation 


[Ref. 


1] : 






[KJ -co 2 [M] 




-[L] 




U 




F 




[*»*] ' 


[K^] 


[0] 






= 


-q 




~p 2 c 2 u) 2 [L] 7 


' [0] r 


[H] -co 2 [M,] 




P 




p 2 C 2 ip 
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where : 



U: 


vector of the nodal values of the components of 
the displacement field, 




vector of the nodal values of the electrial 
potential , 


P: 


vector of the nodal values of the pressure field, 


F: 


vector of the nodal values of the externally 
applied forces, 


T- 


vector of the nodal values of the externally 
applied electrical charges, 


rh 


vector of the nodal values of the integrated 
normal derivative of the externally applied 
pressure on the surface boundary S, 


[* uu ] : 


Finite Element stiffness matrix, 


[*u*] •• 


Finite Element pieozoelectric matrix. 


[K^] : 


Finite Element dielectric matrix, 


[M] : 


Finite Element consistent mass matrix, 


[H] : 


Finite Element fluid pseudo- stiffness matrix, 


[M,] : 


Finite Element consistent fluid pseudo- mass 
matrix, 


[L] : 


Finite Element coupling matrix at the fluid 
structure interface, 


[0] : 


zero matrix, 


w: 


angular frequency, 


P : 


fluid density, 
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fluid sound soeed . 



*: means transposed. 

The results of this analysis for each input frequency 
are the complex displacement, rotation, and electrical 
potential fields at each transducer node, the complex pressure 
field at each fluid node, and the complex electrical impedance 
and admittance. 

2 . Modal Analysis of a Piezoelectric Transducer 

This problem is governed by the equations of motion in 
the elastic and piezoelectric structures, and by Poisson's 
equation in the piezoelectric structures. The matrix equation 
governing this problem is easily obtained from that described 
in the previous section. In a modal analysis there is no fluid 
and there are no external forces applied (the natural boundary 
conditions) , so the third row and column of Equation 13 become 
irrelevant, and F is replaced by 0, resulting in 



[KJ - or [Af] 






~ U " 




" 0 




IK"} 











where the elements are defined by Equation 13. 

In this equation the resonance condition, which 
corresponds to the electrical short-circuit condition, is 
obtained by setting <£=0. The anti - resonance condition, which 
corresponds to the electrical open-circuit condition, is 
obtained by setting g=0 [Ref 15] . 
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The results of this analysis are the eigenfrequencies 
and eigenmoaes. The maximum number of modes, which must be 
specified by the user, is IOC. 

D. STANDARD NETWORK ANALYSIS [Ref. 2] 

The state equations for piezoceramics as stated in 
Equation 9 and 10 can generally be expressed as follows: 

S' = S E t] T - a, E k (15) 

D k = cL T, * e T u E, (16) 

where [T] is the stress tensor, [S] is the strain tensor, E 
is the electric field vector, D is the electric displacement 
vector, [s 5 ] is the constant electric field elastic compliance 
tensor, [d] is the transposed piezoelectric stress tensor, and 
[ e ] is the constant stress permittivity tensor; i, j, k and 
1 can be equal to 1, 2 and 3. The matrix notation has been 

used [Ref . 16] . 

If it is assumed that the ceramic is isotropic in the 
plane transverse to the poling, the piezoelectric equations in 
the cylindrical coordinate system may be written as 
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wnere s L1 «= 


S 33 / 


Sj2'~“ S 


u = ( / u r . 


u 9' 


u z ) , 



. If we let the displacement be expressed as 
then the strain components are 



du 

c - f 

"■ Ir 

_ 3u,_ u, i du r 

V ¥‘T T¥ 



Q 1 du e 

::> M~ — 



&,= 



+ Ur 

r Jd T r 
3 u e i d u. 

~3J + ' 

du. 



"W 



Szz ~~5z 

du. du r 

Srz = ~^¥ + ~Jz 



( 18 ) 



respectively. 
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The equations of motion of the elastic medium in the 
cylindrical coordinate system are 



pu = 


dT rr 






dT._ 


^ _ rn 

1 rr 1 M 


dr 


r* 


dd 


az 
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pu 9 = 






dT„ e 


dT q . 
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3z 
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dT ,Z 
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dd 
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1. Vibration of Unloaded Tube 

For the tangentially polarized ring, the following 
assumptions were made: 

a. the electric field, E is approximated as 

(0, Eq, 0) 

b. the unloaded thin-wall tube vibrates with 
axial symmetry 

c. the stress components, T rz , T zr/ T rr , T rt . 
and T 0r on the outer, inner surfaces and in 
the whole body of the tube are neglected. 

d. T zz =0atz=±l/2 

e . u Q = 0 ; 

f. dT^/de = dT QQ /d 9 = dT zQ /d 6 = 0 
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With the assumptions mentioned, the piezoelectric 
equations may be simplified as 

0 

L > 1 - ~ 'yCijj - ad,, j E „ ) v 20 ) 

[ i r£ ./ (1-J : ) ;(aS W TS„-(d 31 T<7d3 3 ) fij 

- d,,r. ; — € E» 




where Y^ = l/s n “ is the Young's modulus and a = -s,v E: /s 11 £ is the 
Poisson ratio. 



The vibration 
piezoceramic cylindrical 



equations of 
tube can also be 



the thin-walled 
simplified as 



pu r = 




pu 8 = 
pu. = 



0 

dT :z 

dz 
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By combining the piezoelectric equations and the equation 
of motion, the following acoustic wave equation may be 
obtained : 



d 2 u. 




+ K a 2 u. = 0 



where 



r ^2 

CJ 




C"c 




i 2 -i 



( 22 ) 



( 23 ) 
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and lo t = c c /a. c is the acoustic wave speed in che ceramic 
given by: 



c 



C 




(25) 



Applying the boundary conditions, we obtain che following: 

u. = AgSin^z) (2 6) 



Ao = J a E, (27) 



a(d33 + ad 3 i) + f i; (aj)(d3i + gd i 3) 

K4f„( w )*<i 2 ]co S (2v) 



If a and b are the inside and outside radii of the ring 
and N is the number of silvered electrodes, the total electric 
charge on the electrodes is given by 



rl/2 rb 

Q = NQ l = D e drdz 



(29) 
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The electric potential between the two near electrodes is 
V 0 and from Maxwell's eauation, it follows that 



NV\ I 



30) 



ri i i ,t 



Combining the aforementioned relations, the admittance 
equation are obtained as 



r=. 



ju:C 0 k 3] 2 


(f,(a>) (1 + ar) +<j(a + r)y 


z' 

tan ( -i 2C,1 ) 


(cr + r) : 


fj co) 


( 1 -a 2 ) 2 [ (co/ co,) *-l] 


1k,i 


l-a : 






2 J 





tja)C 



(31) 



where 



k 2 



Y E dj\ 



-33 



C. = C 0 [l-(k 31 2 / (l-a : ) ) (l+2crr + r : ) ] 
In (b/a) 



C 0 = N 2 



633 T l ' 



In 

T = ( / dj, ) 



(32) 



These equations are applied in obtaining the results in 
Chapter V, Section E. 
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IV. TWO-DIMENSIONAL FINITE-ELEMENT MODEL 



A. INTRODUCTION 

In order to efficiently model the free flooded ring, as 
shown in Figure 1, the structure may be considered to be an 
axisymmetric body. Appropriate modifications to the material 
properties must be made to model the structure as 
axisymmetric. A continuing problem in the analysis of these 
piezoelectric devices is the determination of these material 
properties . 

A two-dimensional finite element model was developed using 
8 -node isoparametric quadrilateral elements for the 
piezoelectric ceramic, the elastic encapsulant materials, the 
mounting plate, supporting structures, and the fluid. The 
fluid domain extends to lm from the acoustic center. 6 -node 
isoparametric linear elements were used for the solid- fluid 
interface and 3 -node linear elements were used for the dipolar 
radiating damper. These were used to terminate the fluid mesh. 
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B. CHARACTERISTICS OF THE MODEL 



1. Material Properties 

a. Unwrapped Piezoelectric ceramic Ring 

As described before, the transducer has 72 
tangentially poled lead zirconate titanate ceramic staves, 
glued together with epoxy resin, and arranged in a 0.3668 m 
diameter ring. 

In order to simplify the finite element model, an 
equivalent homogeneous piezoelectric ceramic ring with 
material properties determined by a combination of the 
piezoelectric ceramic and epoxy resin properties was used. 
Because the appropriate combination of properties is not 
known, "smeared" material properties for the homogeneous ring 
were obtained by adjusting the circumferential compliance so 
that the resonance frequencies obtained using the ATILA code 
tolerably matched the in- air experimental measurements of the 
resonance frequencies for the unwrapped ring. 

A limitation of an axisymmetric ATILA model is that 
the tangential poling of the segmented ring cannot be modeled 
directly. It has to be modeled as a radially-poled ring as 
illustrated in Figure 9 . 
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MODELED AS RADIALLY-POLED RING 




Figure 9: Modeled as a Radially- Poled Ring 

The transformation from the tangentially to radial 
polarization is accomplished by a suitable exchange of the 
elastic, piezoelectric, and dielectric tensors as described in 
the ATILA manual [Ref. 1] . Since the potential degree of 
freedom is a scalar, the change in polarization direction is 
simply a geometrical transformation. Thus, to model the 
applied potential of "he tangentially poled ring using the 
radially poled ring, the radial potential has to be corrected 
by a factor equal to the ratio of the radial thickness to the 
stave thickness, d/e in Figure 9. Similarly, the results 
computed in a harmonic analysis for the radially-poled 
axisymmetric model have to be modified as follows: (1) the 
displacement and pressure magnitudes must be divided by the 
ratio of the circumferential length of a "smeared" 
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piezoelectric element, which is equal to 1/72 of the ring 
circumferential length, to its thickness, and (2) multiply the 
electrical impedance by the square of the same ratio. This 
ratio is 2.602 for the transducer under study. 

b. Wrapped Piezoelectric Ceramic Ring 

Normally, the fiberglass belt wrapped around the 
ceramic ring (Figure 3) is modeled as thin shell element. 
However, to simplify the model and to reduce the number of 
degrees of freedom, the 0.762 mm thick fiberglass belt in our 
case was not modeled as a separate element. Instead, the 
material properties of the equivalent homogeneous ceramic ring 
were modified. A sample calculation of the approximate new 
material constants is given in Appendix D. These approximate 
values were then used for a modal analysis using ATILA and 
were fine tuned by adjusting first the circumferential 
compliance so that the resonance frequency obtained tolerably 
matched the measured values provided by the manufacturer. Then 
the other compliances were adjusted to obtain a close match of 
the coupling coefficient to the measured value. This scheme is 
considered sufficient for the model even though a scheme 
implemented by McMahon and Armstrong [Ref. 17 ] , which fits the 
resonance and antiresonance frequencies and the electrical 
capacitance taking into consideration the compression bias 
effect of the fiberglass wrapping, might give a more accurate 
result . 
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c. Encapsulant Materials and Mounting Plate 

The encapsulant materials, namely the polyurethane 
and the neoprene, the reflector ring (high carbon steel) , and 
the mounting plate, a fiberglass - epoxy composite, were modeled 
individually as 8 -node quadrilateral elastic elements. The 
material properties were provided by the manufacturer. 

2 . Types of Elements 

The quadratic isoparametric elements listed in Table 
1, which are described in the ATILA User's Manual [Ref. 1] , 
were used: 



Table 1: TYPES OF ELEMENTS USED 



Region 


Element 


Geometry 


Piezoelectric ring 


QUAD08P 


8 -node quadrilateral 


Encapsulant (Boots, 
Spacer, Reflector 
Ring) & Mounting 
Plate 


QUAD08E 


8 -node quadrilateral 


Interface solid- fluid 


LINE06I 


6 -node isoparametric 


Fluid 


QUAD08F 


8 -node quadrilateral 


Radiation surface 


LINE03R 


3 -node linear 
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3. Constraints on Mesh Design [Ref. 1] 

Design of the mesh was guided by the following 
constraints : 

a. Aspect ratio 

The aspect ratio of each element should not be 
greater than 3, although 4 is considered an acceptable, though 
less conservative, value. 

b. Internal angles 

The angles between adjacent sides of quadrilaterals 
should be between 45 and 135 degrees, although 30 degrees and 
150 degrees are considered, respectively, acceptable, though 
less conservative, values. 

c. Element size 

As ATILA utilizes quadratic interpolation 
functions, the size of each element must not be greater than 
one- fourth wavelength at the highest frequency of interest. 

d . In ter el emen t compa tibili ty 

The mesh should be built in such a way that 
adjacent elements have adjoining sides with collocated nodes 
to ensure accurate interpolation at their interfaces. 

e. Radiation boundary elements 

For in- water harmonic analyses (radiation problems) 
the fluid mesh outer limit must be spherical. This is required 
by the radiation damping elements available in the ATILA. code. 
ATILA offers damper elements that are attached to the external 
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surface of the surrounding fluid domain and are designed to 
absorb various components of the radiated multipolar acoustic 
field to overcome the problem of modeling an infinite fluid 
domain. The ATILA offers the monopole and dipole radiation 
damping elements. The latter includes not only the monopole 
term of the radiated field multipolar expansion, but also the 
dipole term. Dipolar damping elements were selected to 
terminate the fluid mesh because, in addition to providing a 
more accurate solution, a smaller mesh can be employed, 
resulting in a lower computational cost. 

A fluid mesh outer limit radius greater than the 
far- field distance is desirable to compute the acoustic source 
pressure level and to compare computed and measured acoustic 
pressure data. The far field criteria for a baffled rigid 
piston-like source is given by [Ref. 15] : 




and 

R ff > L (34) 

where D is the diameter of the ring which is 0.3795 m. 

The minimum wavelength, X, corresponding to the 
maximum frequency of interest is given by: 

X = = 0.375m (35) 
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where c is the speed of sound in water, 1500 m/s, and f is the 
highest frequency of interest. 4000 Hz, giving R ff > 0.384 m. 

The boundary was placed at a radius R equal to 1 m 
from the transducer's acoustical center, which is more than 
2.5 times the far- field limit of the equivalent piston- like 
source at the resonance frequency [Ref. 18] . 
f. Boundary Conditions 

Boundary condition data entries enable the user to 
force, for example, clamped, hinged, or simply- supported 
conditions, or to set the master degrees-of- freedom. These 
conditions are generally associated with symmetry planes or 
axes, electroded or pressure release surfaces. 

In a harmonic analysis of a piezoelectric 
structure, the degree -of- freedom, D, for all the nodes 
associated with the applied potential, V, to the plane or 
line, P, containing the node N has to be a master degree of 
freedom and identical. This is realised by using the following 
line in our data file: 

-N -D -P V 

for example, 

-1 -4 -2 1 

where the degree of freedom 4 indicates electrical potential 
on a plane 2 (normal to the OY axis) passing through the node 
number 1 with an excitation voltage of 1 V. The negative sign 
in front of the 2 indicates that the electrical potential 
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degree of freedom of all nodes on than plane are identical, 
i . e. , IV. 

Since our model is a 2-D model, the degree- of - 
freedom in the plane normal to the OZ axis are all deleted 
using the following data line: 

-N D P 

where D is 3 which corresponds to the displacement in the z- 
direction and P is 3 which corresponds to the plane normal to 
the OZ axis. 

4. Final Mesh Designs 

With the above considerations, an in-air 
axisymmetrical mesh about the OX axis was designed for fitting 
the material properties. Only half of the ring cross-section 
needs to be modeled. The bare ceramic staves without the 
fiberglass wrap were divided into five elements to meet the 
aspect ratio requirement of the ATILA code. The in -air model 
is as shown in Figure 10. 
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Figure 10: In-Air Model 



Using the pre- processor mesh generator MOSAIQUE [Ref. 1] , 
the in- water mesh without the mounting plate as shown in 
Figure 11 was generated. The transducer mesh consists of 40 
solid elements, 21 interface elements, 166 fluid elements, and 
15 radiating elements, totaling 2884 DOF. This mesh was also 
used co investigate effects of encapsulant material changes. 

The final mesh to be developed was the in- water model with 
the mounting plate. The total mesh is shown in Figure 12 and 
an enlarged view of the solid mesh with the mounting plate is 
as shown in Figure 13. There were a total of 87 solid 
elements, 48 interface elements, 406 fluid elements and 35 
radiating elements with a total of 6600 DOF. The input data 
file of this final mesh is given in Appendix A. 
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IN -WATER 
F.E. MODEL 
WITH 

MOUNTING 

PLATE 



Figure 12: In-Water Model With Mounting Plate 
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Figure 13 : The Enlarged View 



V. RESULTS 



A. MATERIAL PROPERTIES OF EQUIVALENT HOMOGENEOUS RING 

The material properties of the equivalent homogeneous 
piezoelectric ceramic ring, unwrapped and wrapped, were 
determined by modal analysis using ATILA. This analysis 
corresponds to a free vibration problem in which the 
eigenf requencies and eigenmodes are computed. 

1. Unwrapped Equivalent Homogeneous Ring 

First, a simple transducer model of only the bare 
ceramic ring without the fiberglass wrap, as shown in Figure 
10, was designed. A modal analysis of this model was conducted 
using book values for the piezoelectric properties. Table 2 
lists the calculated resonance and antiresonance frequencies 
and the coupling coefficient for the first ten modes of 
vibration. The coupling coefficient is calculated using the 
relationship [Ref. 19]: 




where k is the coupling coefficient, f r is the resonance 
frequency and f a is the antiresonance frequency. 
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TABLE 2: RESULTS OF THE MODAL ANALYSIS OF UNWRAPPED 



PIEZOELECTRIC RING USING BOOK VALUES FOR 
THE MATERIAL PROPERTIES 



No. 


RESONANCE 


ANTIRESONANCE 


COUPLING 




( HZ ) 


( HZ ) 


COEFF . , k 




x . 0 6 0 “ J5 


0 . 86e- 3 


- 0 


2 . 


2491 


3287 


0.652 


3 . 


2499 


2499 


0 


4. 


2611 


2611 


0 


5. 


3287 


3488 


0.335 


6. 


4889 


4889 


0 


7. 


7393 


7393 


0 


8. 


10402 


10446 


0.092 


9. 


10728 


10728 


0 


10 . 


14900 


14900 


0 



The ring mode of vibration (mode 2) is the most 
important in operation and is shown in Figure 14 . The dashed 
lines in Figure 14 correspond to the rest position, and the 
deformed shape is in solid lines. 
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X 




To obtain the smeared properties of the unwrapped ring 
assembly, the circumferential compliance s 33 was adjusted so 
that the ring mode resonance frequency obtained for the 
unwrapped ceramic ring in air using the ATILA model matched 
the manufacturer's measured values. 

Table 3 compares the results obtained using typical 
elastic constant values provided in handbooks and the adjusted 
values . 
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TABLE 3 : RESULTS FOR RING MODE OF UNWRAPPED CERAMIC RING 



USING BOOK VALUES AND ADJUSTED VALUES 





MEASURED 


USING BOOK 
VALUES 


ADJUSTED 

VALUES 


Resonance 


2356 Hz 


2491 Hz 


2356 Hz 






(5.7%) 


(0.0%) 


Anti resonance 


4224 Hz 


3287 Hz 


3177 Hz 


Coupling 


0.665 


0 . 652 


0 .671 


Coefficient, k \ 






(0.9%) 



Note: The number in brackets are the percentage error. 

A comparison of the value for the ring mode frequency 
for the ATILA 2-D in-air model for the unwrapped ceramic using 
the adjusted values for elastic constants and the in-air 
experimental results provided by the manufacturer shows the 
expected match of resonance frequency. The computed coupling 
coefficient, however, was 0.9% higher than the measured value. 
The difference is due to the fact that only the elastic 
constant s 33 was adjusted. (Several ATILA modal analysis runs 
were made with adjustments made to the elastic constants s u 
and s 12 just to see the effect on the ring resonance frequency 
and coupling coefficient. The resonance frequency changed only 
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very slightly and the change to the coupling coefficient was 
also very small.; It was felt that a 0.9% difference is 
acceptable for our purposes. For our results, the value of s, 3 
was increased by 12 percent i.e., more compliance; from the 
book value of 0.1539E-10 to 0.1724E-10. This is consistent 
with our expectation that a ring made of glued segments is 
more compliant than a solid piezoelectric ceramic ring. This 
value of s 33 was used for calculating the first -guess smeared 
material properties for the fiberglass wrapped ceramic ring as 
described below. 

2. Wrapped Equivalent Homogeneous Ring 

Smeared material properties for the homogeneous 
equivalent wrapped ceramic ring assembly were obtained by a 
four-step process. First, the coarse values for the elastic 
constants were calculated. A sample calculation is given in 
Appendix D. The circumferential compliance s 33 is then adjusted 
so that the calculated ring resonance frequency matched the 
measured frequency. The other elastic constants were then 
adjusted by a common factor such that the computed coupling 
coefficient matched that of the measured value. The results 
for the ring resonance frequency and coupling coefficient are 
given in Table 4: 
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TABLE 4 : RESULTS OF THE WRAPPED EQUIVALENT CERAMIC RING 





MEASURED 


USING VALUES 


ADJUSTED 






OBTAINED 


VALUES 




1 


EARLIER 












RESONANCE 


2462 Hz 


2442 Hz 


2465 Hz 


FREQUENCY 




(0.8%) 


(0.1%) 


COUPLING 


0 . 646 


0 . 65 


0 . 644 


COEFFICIENT 




(0.62%) 


(0.3%) 



Note: The number in brackets are the percentage error. 

Table 5 lists the values of the elastic constants 
calculated as described in Appendix D, and the final adjusted 
values. These latter values were then used for in-water 
harmonic analyses: 
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TABLE 5: ELASTIC CONSTANTS USED FOR IN-WATER HARMONIC 
ANALYSES 



COMPLIANCE 


COMPUTED AS IN 


ADJUSTED m : /N) 


CONSTANTS 


APPENDIX D \,rtr/N) 




S j j = S 22 


1 . 142e- 11 


1 . 342e- 11 


S j s *> 1 


- 4 . 9e - 12 


-4 . 05e- 12 


S I 3 =S}, 


-5 . 3e - 12 


- 4 . 43e - 12 


S^2 = £>32 


- 5 . 3e - 12 


-4 . 43e- 12 


S 33 


1 . 722e- 11 


1 . 702e- 11 


^ 44=^55 


3 . 9e - 11 


3 . 9e - 11 


S 6 6 


3 . 25e- 11 


3 . 58e- 11 



The other material constants, such as the dielectric 
constants and the piezoelectric constants, were not adjusted. 
The effective density of the homogeneous ring incorporating 
the fiberglass wrap was obtained by dividing the total mass of 
the ceramic and fiberglass by the volume of the homogeneous 
ceramic ring. 



45 



B. IN-WATER HARMONIC ANALYSES 



A harmonic analysis corresponds to a forced vibration 
problem, the excitation being the voltage applied across the 
electrical terminals of the transducer. The applied voltage in 
our case is 1 Vrms . Internal material losses for the 
encapsulant materials were included in this model (but not for 
the ceramic ring assembly) . These were obtained from the loss 
tangents for the respective materials . Values of 0.4 and 0.7 
were used for the loss tangent for polyurethane and neoprene, 
respectively [Ref. 20]. 

1. Mounting Plate Not Included 

Using the material property values obtained from in- 
air modal analyses, described earlier, an in- water model 
without the mounting plate was first analyzed. The mesh is 
shown in Figure 11. 

The plot in Figure 15 depicts the transmitting voltage 
response curve obtained by ATILA for this model along with the 
corresponding manufacturer's data for the actual transducer 
(with the mounting plate) . 

The computed transmitting voltage response was 
generally of the same form except for the lower frequency 
region between 1000 to 1750 Hz. The discrepancy in this region 
is probably due to not including the mounting plate in this 
model. The mounting plate is expected to change the 
transmitting voltage response at the cavity resonance as well 
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as the cavity resonance frequency itself. The TVR at the ring 
resonance frequency was calculated to be about 2 dB lower than 
that measured at the same frequency. 




FREQUENCY 



Figure 15: Comparison Of The Measured Transmitting Voltage 
Response With That Computed For In-water Model 
without Mounting Plate 
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2 . Mounting Plate Included 

An in-water ATILA model including the mounting plate 
was then developed. The mesh are shown in Figures 12 and 13. 
A comparison of the calculated transmitting voltage response 
with the measured response is shown in Figure 16 below: 




FREQUENCY 



Figure 16: Comparison Of Measured TVR With That Calculated 
For In-water Model with Mounting Plate 
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The agreement between the calculated and measured TVR 
is about the same as the previous model. It is greatly 
improved in the cavity mode region and the ring resonance 
frequency calculated is exactly the same as the measured 
values. However, there is a deep dip at 1500 Hz. If this dip 
can be made shallower, then the overall agreement would be 
excellent. It is hoped that the results will be better when 
accurate material property values are available for the 
polyurethane and neoprene encapsulants . The variation in TVR 
with encapsulant material properties is explored in the 
Section D. 

C. COMPARISON OF CAVITY MODE FREQUENCY [Ref. 3] 

The results of the previous section for the cavity mode 
frequency can be compared with the approximate analytical 
formulae given by G.W. McMahon and present previously in 
Chapter III. 

Using relationships given by G.W. McMahon, the first 
cavity mode can be calculated from Equation 2. The speed of 
sound in the water column, c 0 is first computed using Equation 
5 presented in Chapter III. With B equal to 2.18 x 10 9 N/m : and 
I 33 to 5.875 x 10 10 N/irr, c 0 is found to be 846.06 m/s. Q is 
found to be 1.075 and oi is 0.519. With these values, the 
cavity resonance frequency, is 7320.7 rad/s or f c is 1165 
Hz. The relative error between this value and the value 
obtained from the modeling (1250 Hz) is 7%. 
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D. EFFECT OF ENCAPSULANT MATERIAL PROPERTIES CHANGE 

The ATILA model, without the mounting plate, was used to 
explore the effect of changing the material properties of the 
encapsuiant materials, namely the polyurethane around the 
ceramic and the neoprene boots . 

Figures 17 and 18 show the results of the TVR computed for 
different values of the Young"' s modulus for neoprene and 
polyurethane, respectively. The curves are plotted over the 
range of values of Young's modulus quoted in the 
manufacturer's literature. 



TRANSMITTING VOLTAGE RESPONSE 

Varying Polyurethane's Young's Modulus 




FREQUENCY 

Neoprene: Y=1E7 p=0.499 rho=1482 



Figure 17: Effects of Varying Properties of Polyurethane 
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From Figure 17, it can be seen that the best agreement 
with the measured result is for the curve with Y = 4.14 x 10 6 
N/m : and in Figure 18, it was the curve with Y = lx 10 7 N/m : . 
These two values were used for the finite element in-water 
harmonic analyses presented in the previous sections. 

It is observed from the two figures that the value of 
Young's modulus for polyurethane has much greater influence on 
the TVR than that for neoprene. No noticeable change in ring 
resonance frequency was observed in any case. The transmitting 
voltage response was found to increase with a decrease in the 
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Young's modulus of polyurethane. It was, however, the reverse 
m the case of the neoprene. 

E. THE ANALYTICAL MODEL 
1. The Unloaded Case 

The expressions developed based on work by Hong- zhang 
Wang [Ref. 2] were implemented using the MATLAB program. The 
program listing is provided in Appendix B. The effective 
elastic constant values obtained for the unencapsulated, 
fiberglass wrapped ceramic ring assembly of the finite element 
model were used for the analytical model. The admittance 
circle of the unencapsulated fiberglass wrapped ring is 
graphically compared to the measured values in Figure 19 . 
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Figure 19: Comparison of Analytical with Measured Results of 
Wrapped Ring In Air 
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The resonance frequency obtained using the analytical 
method (2420 Hz) was approximately 2 percent lower than the 
measured ^2462 Hz) and the value obtained using the finite 
element code 2465 Hz). 

Wang's approach is quicker and easier to use compared 
to finite element modeling as only the circumferential 
compliance s 33 and the Poisson ratio are required for the 
modeling. It is especially useful when only measured 
admittance circle and handbook material property values are 
available. The expressions developed, however, do not take 
into consideration the encapsulant materials. 

Also, Wang's approach could not be applied to the 
piezoelectric thin ring transducer in water, as it was not 
possible to obtain expressions for the radiation impedances 
presented to the inner and outer surfaces of the ring 
transducer. 
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VI. CONCLUSIONS 



A successful two-dimensional finite element model of a 
free- flooded, tangentially polarized piezoelectric ring sonar 
transducer was built. The model includes internal losses of 
the encapsulant materials. Although the model includes many 
necessary simplifications to handle the problem in the 
available MICROVAX VMS system, they were designed in such a 
way that the transmitting voltage response results obtained 
from the in-water harmonic analyses differs by only 2 dB from 
the measured value at the ring resonance frequency. In the 
region of the cavity mode, the model including the mounting 
plates has better agreement except for a dip at 1500 Hz. It 
was able to accurately predict the ring resonance frequency. 

The model built (without the mounting plate) was used to 
investigate the effect on the acoustical properties of the 
transducer due to a change in the properties of the 
encapsulant materials used in manufacturing. It was observed 
that the Young's modulus of the polyurethane has a significant 
influence on the computed transmitting voltage response. 
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For the analytical model, only the case of the transducer 
in air was developed. The resonance frequency was found to be 
only 2 percent lower than the measured value. 
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VII. SUGGESTIONS FOR FUTURE WORK 



Suggestions for future work include: 

1. Obtain samples of the actual encapsulant materials 
used for testing co determine the actual material 
properties. With these values, the model can then be 
refined. A harmonic analysis should then be performed 
and the transmitting voltage response results obtained 
compared with measured values. 

2. With known values of the material properties, the 
effects of changes m encapsulant materials should 
again be investigated using the in-water model with 
mounting plate. 

3 . Finally, it would be very helpful and interesting 
to obtain an analytical expression for the radiation 
impedance of a free -flooded ring transducer. With this 
and the results of Hong-zhang Wang [Ref. 2], a more 
complete analytical model for the transducer in water 
can be produced. The results obtained with this model 
should then be compared to those obtained using ATILA 
and the measured values . 

4. Study of array element interaction effects should 
be made using these models to develop techniques for 
array performance modeling. 
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APPENDIX A 



INPUT DATA FILE FOR HARMONIC FINITE ELEMENT MODEL 



In-Water Harmonic Analysis of a free-flooded piezoelectric ring 

* sonar transducer 

* Prepared By MAJ Tiong Beng Tay, 14 Feb 1993 

* 

* Navy Type I ceramic * 

* Encapsulant Material: Polyurethane & Neoprene * 

* Tangentially Polarized but modeled as radially-poled * 

* Radius of dipolar damper = 1 m * 



RADIATION DIPOLAR 
ANALYSIS HARMONIC 
SKYLINE COMPLEX 
PRECISION DOUBLE 
CLASS AXISYMME 
LCPDDC = 7 

NLOAD =40 * The number of loading cases must be less than or equal 

* to 100. The minimum is the number of frequencies 

* multiplied by 2 because we are analyzing a radiating 

* transducer 

FREQUENCY 

0 . 225E+04 0 . 250E+04 0.275E+04 0.300E+04 0.325E+04 0.35E4 & 

0.375E4 0.4E4 

GEOMETRY POLARIZA CARTES IA 
1 

0 . 9 00E + 02 0 . 000E+00 0.180E+03 

GEOMETRY 

2 

0.762E-03 0 . 000E+00 

3 

0 .100E+01 



MATERIAL 
PZT4TA 
0 . 000E+00 
0 . 134E-10 
-0 .443E-11 
-0 .405E-11 
0. 000E+00 
0 . 000E+00 
0 . 000E+00 
0 . 000E+00 
0 .000E+00 
-0 . 135E-09 
0 .624E-08 
0 .OOOE+OO 
0 .OOOE+OO 

FGLR 

0 . 587E+11 
P1590 
0 .100E+09 



0 .000E+00 
-0 .443E-11 
0 .170E-10 
-0 .443E-11 
0 . 000E+00 
0. 000E+00 
0 .000E+00 
0 .OOOE+OO 
0 . OOOE+OO 
0 .300E-09 
0. 000E+00 
0 .624E-08 
0 .000E+00 



0.280E+00 
0 .490E+00 



0 . 700E+04 
-0 .405E-11 
-0 .443E-11 
0 .134E-10 
0 .000E+00 
0 .000E+00 
0 .000E+00 
0 .000E+00 
0. 000E+00 
-0 .135E-09 
0 . 000E+00 
0 .000E+00 
0 .570E-08 
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890 


993 


937 


993 


1098 


1041 


1098 


1202 


1145 


1202 


1306 


1249 


1306 


1410 


1353 


1410 


1514 


1457 


1514 


1595 


1562 


582 


691 


636 


691 


795 


740 


795 


900 


842 


900 


1001 


942 


1001 


1103 


1044 


1103 


1208 


1148 


1208 


1312 


1252 


1312 


1416 


1356 


1416 


1521 


1460 


1521 


1601 


1564 


590 


698 


640 


698 


803 


746 


803 


907 


849 


907 


1009 


947 


1009 


1109 


1047 


1109 


1214 


1151 


1214 


1318 


1255 


1318 


1422 


1359 


1422 


1526 


1464 


1526 


1605 


1568 


600 


705 


645 


705 


809 


749 



552 


658 


614 


658 


767 


719 


767 


370 


829 


370 


976 


932 


976 


1081 


1036 


1081 


1186 


1141 


1186 


1290 


1245 


1290 


1394 


1349 


1394 


1498 


1453 


1498 


1581 


1557 


558 


663 


619 


663 


772 


726 


772 


876 


832 


876 


981 


934 


981 


1086 


1039 


1086 


1190 


1143 


1190 


1294 


1247 


1294 


1398 


1351 


1398 


1502 


1455 


1502 


1585 


1559 


564 


672 


627 


672 


776 


734 


776 


882 


836 


882 


986 


937 


986 


1090 


1041 


1090 


1194 


1145 


1194 


1298 


1249 


1298 


1402 


1353 


1402 


1506 


1457 


1506 


1589 


1562 


572 


679 


636 


679 


785 


740 


785 


888 


842 


888 


990 


942 


990 


1094 


1044 


1094 


1198 


1148 


1198 


1303 


1252 


1303 


1408 


1356 


1408 


1512 


1460 


1512 


1593 


1564 


580 


687 


640 


687 


792 


746 


792 


896 


849 


896 


998 


947 


998 


1101 


1047 


1101 


1205 


1151 


1205 


1309 


1255 


1309 


1413 


1359 


1413 


1517 


1464 


1517 


1598 


1568 


585 


695 


645 


695 


800 


749 


800 


903 


854 


903 


1004 


956 


1004 


1107 


1054 


1107 


1211 


1155 


1211 


1316 


1259 


1316 


1420 


1363 


1420 


1524 


1467 


1524 


1603 


1571 


596 


701 


651 


701 


807 


755 



90 



- 320 * 


803 


907 


* 321 * 


907 


1009 


* 322 * 


1009 


1109 


* 323 * 


1109 


1214 


* 3 24 * 


1214 


1318 


- 325 * 


1318 


1422 


* 326 * 


1422 


1526 


* 327 * 


1526 


1605 


* 328 * 


600 


705 


* 329 * 


705 


809 


* 330 * 


809 


915 


* 331 * 


915 


1017 


* 332 * 


1017 


1118 


* 333 * 


1118 


1220 


* 334 * 


1220 


1324 


* 335 * 


1324 


1428 


* 336 * 


1428 


1532 


* 337 * 


1532 


1609 


* 338 * 


625 


724 


* 339 * 


724 


822 


* 340 * 


822 


925 


* 341 * 


925 


1029 


* 342 * 


1029 


1131 


* 343 * 


1131 


1229 


* 344 * 


1229 


1333 


* 345 * 


1333 


1437 


* 346 * 


1437 


1541 


* 347 * 


1541 


1614 


* 348 * 


608 


710 


* 349 * 


710 


813 


* 350 * 


813 


916 


* 351 * 


916 


1021 


* 352 * 


1021 


1120 


* 353 * 


1120 


1222 


* 354 * 


1222 


1326 


* 355 * 


1326 


1430 


* 356 * 


1430 


1534 


* 357 * 


1534 


1618 


* 358 * 


603 


703 


* 359 * 


703 


806 


* 360 * 


806 


910 


* 361 * 


910 


1011 


* 362 * 


1011 


1110 


* 363 * 


1110 


1213 


* 364 * 


1213 


1317 


* 365 * 


1317 


1421 


* 366 * 


1421 


1525 


* 367 * 


1525 


1622 


* 368 * 


594 


693 


* 369 * 


693 


793 


* 370 * 


793 


892 


* 371 * 


892 


992 


* 372 * 


992 


1096 


* 373 * 


1096 


1200 


* 374 * 


1200 


1302 


* 375 * 


1302 


1406 


* 376 * 


1406 


1510 


* 377 * 


1510 


1626 


* 378 * 


583 


677 


* 379 * 


677 


769 


* 380 * 


769 


868 


* 381 * 


868 


972 



309 


915 


854 


915 


1017 


956 


1017 


1118 


1054 


1118 


1220 


1155 


1220 


1324 


1259 


1324 


1428 


1363 


1428 


1532 


1467 


1532 


1609 


1571 


625 


724 


651 


724 


822 


755 


822 


925 


858 


925 


1029 


961 


1029 


1131 


1060 


1131 


1229 


1157 


1229 


1333 


1261 


1333 


1437 


1365 


1437 


1541 


1469 


1541 


1614 


1573 


608 


710 


670 


710 


813 


764 


813 


916 


866 


916 


1021 


969 


1021 


1120 


1071 


1120 


1222 


1167 


1222 


1326 


1266 


1326 


1430 


1370 


1430 


1534 


1474 


1534 


1618 


1578 


603 


703 


654 


703 


806 


757 


806 


910 


861 


910 


1011 


963 


1011 


1110 


1064 


1110 


1213 


1160 


1213 


1317 


1263 


1317 


1421 


1367 


1421 


1525 


1471 


1525 


1622 


1575 


594 


693 


650 


693 


793 


752 


793 


892 


855 


892 


992 


957 


992 


1096 


1055 


1096 


1200 


1154 


1200 


1302 


1258 


1302 


1406 


1362 


1406 


1510 


1466 


1510 


1626 


1570 


583 


677 


643 


677 


769 


744 


769 


868 


845 


868 


972 


944 


972 


1077 


1043 


1077 


1181 


1147 


1181 


1285 


1251 


1285 


1389 


1355 


1389 


1493 


1459 


1493 


1630 


1561 


573 


655 


633 


655 


756 


728 


756 


860 


823 


860 


966 


929 



807 


911 


858 


911 


1013 


961 


1013 


1114 


1060 


1114 


1217 


1157 


1217 


1321 


1261 


1321 


1425 


1365 


1425 


1529 


1469 


1529 


1607 


1573 


611 


713 


670 


713 


816 


764 


816 


920 


866 


920 


1023 


969 


1023 


1124 


1071 


1124 


1224 


1167 


1224 


1328 


1266 


1328 


1432 


1370 


1432 


1536 


1474 


1536 


1611 


1578 


616 


716 


654 


716 


817 


757 


817 


922 


861 


922 


1025 


963 


1025 


1127 


1064 


1127 


1226 


1160 


1226 


1330 


1263 


1330 


1434 


1367 


1434 


1538 


1471 


1538 


1616 


1575 


604 


706 


650 


706 


811 


752 


811 


914 


855 


914 


1016 


957 


1016 


1117 


1055 


1117 


1219 


1154 


1219 


1323 


1258 


1323 


1427 


1362 


1427 


1531 


1466 


1531 


1620 


1570 


599 


699 


643 


699 


801 


744 


801 


902 


845 


902 


1002 


944 


1002 


1104 


1043 


1104 


1207 


1147 


1207 


1311 


1251 


1311 


1415 


1355 


1415 


1518 


1459 


1518 


1624 


1561 


588 


685 


633 


685 


783 


728 


783 


881 


823 


881 


980 


929 


980 


1079 


1034 


1079 


1183 


1138 


1183 


1287 


1242 


1287 


1391 


1346 


1391 


1495 


1450 


1495 


1628 


1554 


578 


665 


622 


665 


759 


711 


759 


863 


812 


863 


970 


919 



91 



* 382 * 


972 


1077 


966 


* 383 * 


1077 


1181 


1072 


- 384 * 


1181 


1285 


1177 


* 385 * 


1285 


1389 


1281 


* 386 * 


1389 


1493 


1385 


* 387 * 


1493 


1630 


1489 


* 388 * 


573 


655 


541 


* 389 * 


655 


756 


641 


* 390 * 


7 56 


860 


747 


* 391 * 


860 


966 


850 


* 392 * 


966 


1072 


958 


* 393 * 


1072 


1177 


1066 


* 394 * 


1177 


1281 


1173 


* 395 * 


1281 


1385 


1278 


* 396 * 


1385 


1489 


1381 


* 397 * 


1489 


1634 


1485 


* 398 * 


541 


641 


510 


* 399 * 


641 


747 


624 


* 400 * 


747 


850 


736 


* 401 * 


850 


958 


841 


* 402 * 


958 


1066 


949 


* 403 * 


1066 


1173 


1056 


* 404 * 


1173 


1278 


1166 


* 405 * 


1278 


1381 


1272 


* 406 * 


1381 


1485 


1376 


* 407 * 


1485 


1637 


1480 


* 408 * 


510 


624 


464 


*4 09 * 


624 


736 


598 


* 410 * 


736 


841 


715 


* 411 * 


841 


949 


831 


* 412 * 


949 


1056 


941 


* 413 * 


1056 


1166 


1051 


* 414 * 


1166 


1272 


1162 


* 415 * 


1272 


1376 


1270 


* 416 * 


1376 


1480 


1374 


* 417 * 


1480 


1640 


1478 


* 418 * 


464 


598 


448 


* 419 * 


598 


715 


591 


* 420 * 


715 


831 


708 


* 421 * 


831 


941 


824 


* 422 * 


941 


1051 


939 


* 423 * 


1051 


1162 


1048 


* 424 * 


1162 


1270 


1159 


* 425 * 


1270 


1374 


1268 


* 426 * 


1374 


1478 


1372 


* 427 * 


1478 


1643 


1476 


* 428 * 


404 


414 


448 


* 429 * 


414 


469 


464 


* 430 * 


383 


400 


404 


* 431 * 


400 


447 


414 


* 432 * 


360 


379 


383 


* 433 * 


379 


433 


400 


* 434 * 


352 


369 


360 


* 435 * 


369 


418 


379 


* 436 * 


351 


368 


335 


* 437 * 


368 


417 


348 


* 438 * 


335 


348 


322 


* 439 * 


348 


398 


337 


* 440 * 


322 


337 


315 


* 441 * 


337 


381 


329 


* 442 * 


315 


329 


311 


* 443 * 


329 


371 


323 



1034 


970 


1075 


1026 


1138 


1075 


1179 


1136 


1242 


1179 


1283 


1240 


1346 


1283 


1387 


1344 


1450 


1387 


1491 


1448 


1554 


1491 


1632 


1552 


622 


559 


646 


587 


711 


646 


750 


692 


812 


750 


856 


799 


919 


856 


962 


906 


1026 


962 


1068 


1018 


1136 


1068 


1176 


1129 


1240 


1176 


1280 


1238 


1344 


1280 


1383 


1342 


1448 


1383 


1487 


1446 


1552 


1487 


1635 


1550 


587 


523 


632 


567 


692 


632 


739 


674 


799 


739 


846 


784 


906 


846 


954 


895 


1018 


954 


1062 


1007 


1129 


1062 


1170 


1121 


1238 


1170 


1275 


1235 


1342 


1275 


1379 


1339 


1446 


1379 


1482 


1443 


1550 


1482 


1638 


1547 


567 


483 


606 


539 


674 


606 


721 


653 


784 


721 


835 


761 


895 


835 


945 


884 


1007 


945 


1053 


999 


1121 


1053 


1163 


1115 


1235 


1163 


1271 


1233 


1339 


1271 


1375 


1338 


1443 


1375 


1479 


1442 


1547 


1479 


1641 


1546 


539 


452 


593 


532 


653 


593 


709 


647 


761 


709 


825 


760 


884 


825 


940 


878 


999 


940 


1049 


995 


1115 


1049 


1161 


1112 


1233 


1161 


1269 


1230 


1338 


1269 


1373 


1337 


1442 


1373 


1477 


1441 


1546 


1477 


1645 


1545 


408 


427 


441 


452 


439 


441 


489 


483 


385 


393 


412 


408 


416 


412 


460 


439 


365 


373 


389 


385 


402 


389 


437 


416 


354 


358 


375 


365 


391 


375 


429 


402 


353 


344 


356 


339 


390 


356 


410 


362 


339 


331 


341 


325 


362 


341 


387 


350 


325 


318 


333 


320 


350 


333 


377 


346 


320 


313 


327 


316 


346 


327 


367 


342 



1072 

1177 

1281 

1385 

1489 

1634 

641 

747 

850 

958 

1066 

1173 

1278 

1381 

1485 

1637 

624 

736 

841 

949 

1056 

1166 

1272 

1376 

1480 

1640 

598 

715 

831 

941 

1051 

1162 

1270 

1374 

1478 

1643 

591 

708 

824 

939 

1048 

1159 

1268 

1372 

1476 

1648 

464 

510 

414 

469 

400 

447 

379 

433 

348 

398 

337 

381 

329 

371 

323 

363 



92 



*44 4 * 


311 


323 


314 


*44 5 * 


323 


363 


328 


*44 6 * 


314 


328 


321 


*44 7 * 


328 


370 


336 


*44 8 * 


321 


336 


334 


* 449 * 


336 


380 


347 


* 450 * 


334 


347 


352 


* 451 * 


347 


397 


369 


* 452 * 


359 


378 


351 


* 453 * 


378 


432 


368 


* 454 * 


382 


399 


359 


* 455 * 


399 


446 


378 


* 456 * 


403 


413 


382 


* 457 * 


413 


468 


399 


* 458 * 


449 


465 


403 


* 459 * 


465 


511 


413 


* 460 * 


474 


490 


449 


* 461 * 


490 


527 


465 


* 462 * 


1613 


1544 


1617 


* 463 * 


1544 


1440 


1539 


* 464 * 


1440 


1336 


1435 


* 465 * 


1336 


1234 


1331 


* 466 * 


1234 


1135 


1227 


* 467 * 


1135 


1032 


1128 


* 468 * 


1032 


930 


1027 


* 469 * 


930 


834 


923 


* 470 * 


834 


737 


819 


* 471 * 


737 


634 


722 


* 472 * 


1617 


1539 


1621 


* 473 * 


1539 


1435 


1533 


* 474 * 


1435 


1331 


1429 


* 475 * 


1331 


1227 


1325 


* 476 * 


1227 


1128 


1221 


*477 * 


1128 


1027 


1119 


* 478 * 


1027 


923 


1020 


* 479 * 


923 


819 


917 


* 480 * 


819 


722 


814 


* 481 * 


722 


623 


714 


* 482 * 


1621 


1533 


1625 


* 483 * 


1533 


1429 


1520 


* 484 * 


1429 


1325 


1417 


* 485 * 


1325 


1221 


1313 


* 486 * 


1221 


1119 


1210 


* 487 * 


1119 


1020 


1106 


* 488 * 


1020 


917 


1005 


* 489 * 


917 


814 


905 


* 490 * 


814 


714 


804 


* 491 * 


714 


617 


702 


* 492 * 


1625 


1520 


1629 


* 493 * 


1520 


1417 


1496 


* 494 * 


1417 


1313 


1392 


* 495 * 


1313 


1210 


1288 


* 496 * 


1210 


1106 


1184 


* 497 * 


1106 


1005 


1083 


* 498 * 


1005 


905 


983 


* 499 * 


905 


804 


885 


* 500 * 


804 


702 


788 


* 501 * 


702 


605 


689 


* 502 * 


1629 


1496 


1633 


* 503 * 


1496 


1392 


1492 


* 504 * 


1392 


1288 


1388 


* 505 * 


1288 


1184 


1284 



316 


312 


326 


319 


342 


326 


366 


345 


319 


317 


332 


324 


345 


332 


376 


349 


324 


330 


340 


338 


349 


340 


386 


361 


338 


343 


355 


354 


361 


355 


409 


391 


364 


357 


374 


353 


401 


374 


428 


390 


384 


372 


388 


364 


415 


388 


436 


401 


407 


392 


411 


384 


438 


411 


459 


415 


453 


426 


440 


407 


484 


440 


488 


438 


478 


463 


477 


453 


509 


477 


518 


484 


1579 


1615 


1543 


1576 


1475 


1543 


1439 


1472 


1371 


1439 


1335 


1368 


1267 


1335 


1232 


1264 


1174 


1232 


1134 


1164 


1074 


1134 


1031 


1069 


973 


1031 


928 


967 


875 


928 


828 


864 


782 


828 


731 


762 


681 


731 


630 


668 


1576 


1619 


1537 


1574 


1472 


1537 


1433 


1470 


1368 


1433 


1329 


1366 


1264 


1329 


1225 


1262 


1164 


1225 


1126 


1158 


1069 


1126 


1024 


1063 


967 


1024 


921 


964 


864 


921 


818 


862 


762 


818 


717 


758 


668 


717 


621 


656 


1574 


1623 


1528 


1567 


1470 


1528 


1424 


1463 


1366 


1424 


1320 


1360 


1262 


1320 


1216 


1256 


1158 


1216 


1113 


1152 


1063 


1113 


1014 


1050 


964 


1014 


913 


953 


862 


913 


810 


853 


758 


810 


707 


751 


656 


707 


613 


652 


1567 


1627 


1511 


1555 


1463 


1511 


1407 


1451 


1360 


1407 


1305 


1347 


1256 


1305 


1201 


1243 


1152 


1201 


1097 


1139 


1050 


1097 


996 


1038 


953 


996 


897 


936 


853 


897 


796 


838 


751 


796 


697 


742 


652 


697 


602 


642 


1555 


1631 


1494 


1553 


1451 


1494 


1390 


1449 


1347 


1390 


1286 


1345 


1243 


1286 


1182 


1241 



328 

370 

336 

380 

347 

397 

369 

418 

368 

417 

378 

432 

399 

446 

413 

468 

465 

511 

1539 

1435 

1331 

1227 

1128 

1027 

923 

819 

722 

623 

1533 

1429 

1325 

1221 

1119 

1020 

917 

814 

714 
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1520 

1417 

1313 

1210 

1106 

1005 

905 
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702 

605 

1496 

1392 

1288 

1184 

1083 

983 

885 

788 

689 

597 

1492 

1388 

1284 

1180 



93 



1184 

1083 

983 

385 

788 

689 

1633 

1492 

1388 

1284 

1180 

1076 

971 

367 

766 

678 

1639 

1488 

1384 

1279 

1175 

1067 

959 

352 

745 

638 

1644 

1484 

1380 

1276 

1171 

1061 
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844 

733 

618 

3R 

1648 

1643 

1640 

1637 

1634 

1630 

1626 

1622 

1618 

1614 

1609 

1605 

1601 

1595 

1591 

1587 

1583 

1580 

1584 

1588 

1592 

1596 

1599 

1602 



1083 

983 

385 

788 

689 

597 

1492 

1388 

1284 

1180 

1076 

971 

867 

766 

678 

586 

1488 

1384 

1279 

1175 

1067 

959 

852 

745 

638 

527 

1484 

1380 

1276 

1171 

1061 

951 

844 

733 

618 

490 



1180 

1076 

971 

867 

^66 

678 

1639 

1488 

1384 

1279 

1175 

1067 

959 

852 

745 

638 

1644 

1484 

1380 

1276 

1171 

1061 
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844 

733 

618 

1647 

1481 

1377 

1273 

1168 

1057 

948 

839 

725 

607 



1076 

971 

867 

^66 

678 

586 

1488 

1384 

1279 

1175 

1067 
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638 

527 

1484 

1380 

1276 

1171 

1061 

951 

844 

733 

618 

490 

1481 

1377 

1273 

1168 

1057 

948 

839 

725 

607 
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1139 
1038 
936 
838 
742 
64 2 
1553 
1449 
1345 
1241 
1137 
1033 
926 
821 
729 
635 
1551 
1447 
1343 
1239 
1132 
1019 
908 
798 
688 
577 
1549 
1445 
1341 
1237 
1125 
1012 
899 
787 
671 
556 



1182 

1078 

974 

872 

^78 

682 

1636 

1490 

1386 

1282 

1178 

1073 

965 

359 

753 

649 

1642 

1486 

1382 

1277 

1172 

1065 

955 

348 

738 

629 

1646 

1483 

1378 

1274 

1169 

1059 

950 

840 

730 

609 



1078 

974 

372 

778 

682 

592 

1490 

1386 

1282 

1178 

1073 

965 

859 

753 

649 

555 

1486 

1382 

1277 

1172 

1065 

955 

348 

738 

629 

509 

1483 

1378 

1274 

1169 

1059 

950 

840 

730 

609 

478 



1137 

103*3 

926 

321 

729 

635 

1551 

1447 

1343 

1239 

1132 

1019 

908 

798 

688 

577 

1549 

1445 

1341 

1237 

1125 

1012 

899 

787 

671 

556 

1548 

1444 

1340 

1236 

1123 

1008 

893 

779 

660 

547 



FLUIDE 3 

1643 1645 
1640 1641 
1637 1638 
1634 1635 
1630 1632 
1626 1628 
1622 1624 
1618 1620 
1614 1616 
1609 1611 

1605 1607 

1601 1603 

1595 1598 

1591 1593 

1587 1589 

1583 1585 
1580 1581 

1584 1582 

1588 1586 

1592 1590 

1596 1594 
1599 1597 

1602 1600 

1606 1604 
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*566* 


1606 


*5 6 7* 


1610 


*568* 


1613 


-569* 


1617 


*570* 


1621 


*571* 


1625 


*572* 


1629 


*573* 


1633 


*574* 


1636 


*575 * 


1639 


*576* 


1644 



1610 1608 
1613 1612 
1617 1615 
1621 1619 
1625 1623 
1629 1627 
1633 1631 
1636 1649 
1639 1650 
1644 1642 
1647 1646 



END 



-1 

-3 



-4 



3 

2 S 
-2 1 



* Non existing degree of freedom [2-D) 

* Electrode at 0 V 

* Electrode at 1 V 
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APPENDIX B 



MATLAB PROGRAM FOR ANALYTICAL MODEL OF UNLOADED TRANSDUCER 



% This MATLAB program was written by MAJ Tay Tiong Beng 
% on 1 Feb 1993 to implement the formula developed 
% by Hong -zhang Wang [Ref. 2] for unloaded tangentially 
% polarized ring transducer 



% Initialization 

clear; 
clg ; 

!del res 
!del resl.met 
!del res2.met 
! del res 3 . met 
diary res 

% Material Constants 

j =sqrt ( -1) ; 
rhoc=7002 ; 
sll=0.1702e-10* (1-0 

sl2=-0 . 405e - 11 ; 
nu=- ( sl2 ) / (sll) 
s66=2* (sll - sl2 ) 
e33=0 . 57e- 8 
d31=- 0 . 135e- 9 
d33=0 .3e-9 
Yc=l/sll ; 

F=3 00 0 ; 

% Fluid Constants 

cair=343 ; 
cwater=1500 ; 
ccer=sqrt (Yc/rhoc) ; 

% Dimensions 

a= ( 14 . 44 /2 } *0 . 0254 % Inside radius 

b= (14 . 94/2) *0 . 0254 % outside radius of ceramic stave 

1 = 6 .8*0 .0254 % height 

% Capacitance and Admittance 

N=72; 

Co=N"2*e33*l/ (2*pi) *log(b/a) 
k31s=Yc* (d3l"2) /e33 
t= (d33/d31 ) 

CssCo* (1- (k31s/ (l-nu*2) ) * (l+2*nu*t+t^2 ) ) 



. 0 1 5 * j ) ; % Losses is 1.5% of the real part. The real part 

% is the published value of the elastic constant 

% Poisson Ratio 
% Frequency 
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% Admittance Circle 



j =sqrt ( -1) ; 

wr=ccer/a % Radial Res 

f r=wr / ( 2 *pi ) 
n = 0 ; 

for f =1 : 5 : F 
n=n+l ; 
w=2 *pi ; 
ww=w/wr ; 

WW (1 , n) =ww ; 
j wcs = j w w*Cs ; 

f a= { 1 -nu^2 ) * (w/wr ) ^2 - 1; 

A= j ^w^Co^kSls/ f a ; 

ka2= (w/ccer) A 2 *f a/ ( (w/wr) A 2 - 1) ; 

KA2 (1 ,n) =ka2 ; 
ka=sqrt (ka2 ) ; 

KA ( 1 , n) =ka ; 

B= (f a* (l+nu*t ) +nu* (nu+t) )^2/((l-nu A 2) A 2*( (w/wr) *2 - 1) ) ; 

D = (nu+t) *2/ (1 -nu A 2) ; 

T=A* (B* (tan (0 . 5*ka*l) / (0 . 5*ka*l) ) -D) + j wcs ; 

TT=T/ ( j *wr^Co) ; 

2=1/ TT ; 

Z ( 1 , n) =z ; 

Y ( 1 , n) =T ; 

YY(l,n) =TT ; 

end; 

q=abs (Y) 

G=real (Y) 

S=imag (Y) 

f = 1 ; 5 : F ; 

axis ( ' square' ) ; 

plot (G, S) ; grid; 

title ('ADMITTANCE CIRCLE OF A TANGENTIALLY POLARIZED UNLOADED TRANSDUCER') 
xlabel (' Conductance , G mhos'); 
ylabel ( ' Susceptance , S mhos'); 
meta resl ; 

plot (f ' , G) ; grid; 

title ('CONDUCTANCE AS A FUNCTION OF FREQUENCY'); 

xlabel (' Frequency, Hz'); ylabel (' Conductance , G mhos'); pause; 
meta res2 ; 

plot (f ' , S) ; grid; 

title (' SUSCEPTANCE AS A FUNCTION OF FREQUENCY') ; 

xlabel (' Frequency , Hz'); ylabel (' Susceptance , S mhos'); pause; 

meta res3 ; 

diary off 
end ; 
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APPENDIX C 



TYPICAL VALDES FOR MATERIALS CONSTANTS 



TABLE 1: PHYSICAL DIMENSIONS AND MATERIAL CONSTANTS 



The Transducer 



Inside Radius 


0.1651 m 


Outside Radius 


0.204724 m 


Height 


0.2921 m 


Motor Element: Naw Type 1 Ceramic (Hexaqonal Clas 


Density 


7550 kg/m 3 


Inside Radius 


0.183388 m 


Outside Radius 


0.189738 m 


Elastic Constants (nr/N) 


Piezoelectric Constants 




a tensor (rrr/c) 


s 33 e = 1.539E-11 




s I3 e = -5.3 IE- 12 


g 33 = 26. IE - 3 


S]2 e = -4.05E-12 


g 31 = -11.7E-3 


s u E = 1.22E-11 


g 15 = 40.5E-3 


S44 E = 3.90E-11 




s 66 e = 3.25E-11 


d tensor (m/v) 


Dielectric Constants 


d 33 = 300E-12 




d 3J = -135E-12 


e 33 T = 1.149E-8 


d l5 = 525E-12 


e,, T = 1 . 2963E- 8 




e 33 s = 0.57E-8 




6 U S = 0.624E-8 
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TABLE 1: CONTINUED 



Fiberglass Wrap 



Thickness 

Density 

Longitudinal Modulus 
Transverse Modulus 
Poisson Ratio 



0.762 mm 
254 0 kg/m 3 

8.5E6 psi (5.86E10 Pa) 
2.9E6 psi (2.00E10 Pa) 
0.28 



Encaosulant Material 



Polyurethane 

Density : 0.039 pci ( >=1100 kg/m J ) 
Young's Modulus : 600 psi (4.14 MPa) 
Poisson Ratio : 0.49 



Neoprene 

Density : 0.045 pci ( «1250 kg/m 3 ) 

Young's Modulus : 100-3000 psi 

(0.69-20.69 MPA) 

Poisson Ratio : 0.49 

Reflector Rina : HRS Carbon Steel 

Density : 7700 Kg/m 3 

Young's Modulus : 19.5E10 Pa 
Poisson Ratio : 0.28 



Mounting Plate : Fiberglass -epoxy composite 

Density : 1700 Kg/m 3 

Young Modulus : 2000 psi (13.79 MPa) 

Poisson Ratio : 0.34 
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APPENDIX D 



SAMPLE CALCULATIONS FOR THE EFFECTIVE 
ELASTIC CONSTANTS OF THE FIBERGLASS WRAPPED CERAMIC 



The effective elastic constants of the ring are estimated as 
follows. Consider a section of the fiberglass wrapped ring as two 
rectangular volume elements cemented together as shown in Figure 1 
below. The purpose is to calculate the elastic constants of an 
equivalent ceramic material which will give the same elastic 
properties as the fiberglass wrapped ceramic ring. The equivalent 
ceramic material has the same length, 1, and height, h, as the 
fiberglass and ceramic, and a width equal to the sum of the widths 
of the fiberglass and ceramic. 



RING 





Figure 1 
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The thickness of the fiberglass is t and that of the ceramic is w, 
where t and w are 0.000762m and 0.0064m, respectively. 

Consider first an applied force, F } , in the circumferential 
direction only, as shown. Assuming that the resultant stress, T } , 
is uniform then the strain, S , , in each material is given by: 

S = = s T (1) 



For the ceramic, 



(Ah) c - h c ( s 13 ) c T 3 



( 2 ) 



And for the fiberglass, 

(Ah) g = h g (s l3 ) g T 3 (3) 

There is no reason to assume (Ah) g. equals to (Ah) c . We 
therefore define (Ah ) eff for the equivalent material as the average 
Ah across the composite material, neglecting boundary effects: 

(Ah) c w+ (Ah) g t = (Ah) eff (w+t) (4) 

Then the effective strain in the axial direction, i.e., in the "1" 
direction is given by: 

<S,> , = (Ah) « - . < S .3>.“"< S .3>A j, (5) 

h w+t 

after taking h c = h g = h. 
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is given by: 



Thus the effective compliance, (s ]3 ) e£f , 



(s 



13 1 ejff 



(s l3 ) c w+(s l3 ) ^ t 
w+ c 



( 6 ) 



By considering forces in the other directions, in turn, the 
same result is obtained, with the appropriate identification of 
indices and element dimensions. For example, to compute s i; we will 
consider an applied force in the "2" direction, i.e., the radial 
direction, and calculate the strain in the ”1", or axial, 
direction . 



Using values listed in Appendix C for the fiberglass and 
ceramic, the following were obtained: 



CONSTANTS 


COMPUTED VALUES (m : /N) 


Sjj — £>*> 2 


1.142 x 10' 11 


s n = s 2l 


-4.9 x 10' 12 


S 13 = S 3I 


-5.3 x 10' 12 


S 23 = S 32 


-5.3 x 10' 12 


S 33 


1.702 x lO' 11 


S 44 = S 55 


3.9 x 10'" 


S 66 


3.58 x 10’“ 
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